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The relation of multihadron production in nucleus-nucleus and (anti)proton-proton collisions is 
studied by exploring the collision-energy and centrality dependencies of the charged particle mean 
multiplicity in the measurements to date. The study is performed in the framework of the recently 
proposed effective-energy approach which, under the proper scaling of the collision energy, combines 
the constituent quark picture with Landau relativistic hydrodynamics counting for the centrality- 
defined effective energy of participants. Within this approach, the multiplicity energy dependence 
and the pseudorapidity spectra from the most central nuclear collisions are well reproduced. The 
study of the multiplicity centrality dependence reveals a new scaling between the measured pseu¬ 
dorapidity spectra and the calculations. By means of this scaling, referred to as energy-balanced 
limiting fragmentation scaling, one reproduces the pseudorapidity spectra for all centralities. The 
scaling elucidates some differences in the multiplicity and midrapidity density centrality dependence 
obtained at RHIC and LHC. These findings reveal an inherent similarity in the multiplicity energy 
dependence from the most central collisions and centrality data. Predictions are made for the mean 
multiplicities to be measured in proton-proton and heavy-ion collisions at the LHC. 

PACS numbers: 25.75.Dw, 25.75.Ag, 24.85.-l-p, 13.85.Ni 


I. INTRODUCTION 


Study of global observables of multiparticle production 
and their universality in different types of high-energy 
collisions is of a crucial importance for understanding the 
underlying dynamics of strong interactions. Recently, 
the universality of multiparticle production in nucleus- 
nucleus and hadron-hadron collisions has been reported 
exploiting concept of centrality-defined effective energy 
[l[ employed for the data interpreted in terms of the ap¬ 
proach of the dissipating energy of quark participants 
[3, 0. This approach combines the constituent quark 
picture together with Landau relativistic hydrodynam¬ 
ics and interrelates multihadron production in different 
types of collisions. Within this picture one can success¬ 
fully explain ii the scaling between the charged parti¬ 
cle mean multiplicity in e+e“ and pp/pp collisions and 
the universality of both the multiplicity and the midra¬ 
pidity pseudorapidity density measured in the most cen¬ 
tral nuclear collisions and in e+e“ annihilation The 
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universality of the multihadron production irrespective 
of the collision species, an intrinsic feature of the dis¬ 
rating energy approach, is widely discussed nowadays 

m 

In this paper, in the framework of the approach of the 
dissipating effective energy of constituent quark partici¬ 
pants, or, for brevity, the participant dissipating energy 
(PDE) approach, we extend the previous studies of the 
charged particle mean multiplicity ii to LHC energies. 
We show that the multiplicity energy dependence of 
head-on collisions is well described within the proposed 
approach. In addition, here we study the dependence of 
the multiplicity on the number of (nucleon) participants 
at RHIC and LHC. We introduce a new scaling, referred 
to as the energy-balanced limiting fragmentation scaling, 
which allows us to describe the pseudorapidity density 
spectra independently of the centrality of collisions. 
Using this scaling, a complementarity between the 
multiplicities measured in head-on nuclear collisions and 
obtained from the centrality data is found. The study 
clarifies some differences of the centrality dependence 
of multiplicities measured at RHIC and LHC. Finally, 
predictions are made for the charged particle mean 
multiplicities in pp and heavy-ion collisions at the LHC. 

II. THE PARTICIPANT DISSIPATING ENERGY 
APPROACH 

In this section, we briefly describe the PDE approach, 
as it is proposed in [in. This approach quantifies the 
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process of particle production in terms of the amount of 
energy deposited by interacting constituent quark partic¬ 
ipants inside the small Lorentz-contracted volume formed 
at the early stage of a collision. The whole process of a 
collision is then represented as the expansion of an initial 
state and the subsequent breakup into particles. This ap¬ 
proach resembles the Landau phenomenological hydrody¬ 
namic approach of multiparticle production in relativis¬ 
tic particle interactions [^, which was found to be in 
a good agreement with the multiplicity data in particle 
and nuclear collisions in the wide energy range |lC| . In 
the picture considered here, the Landau hydrodynamics 
is employed in the framework of constituent (or dressed) 

a arks, in accordance with the additive quark model [n^ 
; for recent comprehensive review on soft hadron inter¬ 
actions in the additive quark model, see [l5l | . This means 
the secondary particle production is basically driven by 
the amount of the initial effective energy deposited by 
constituent quarks into the Lorentz-contracted region. In 
pp/pp collisions, a single constituent quark from each nu¬ 
cleon is considered to take part in a collision, and the re¬ 
maining quarks are treated as spectators. The spectator 
quarks do not participate in the secondary particle pro¬ 
duction, but they result in a formation of leading par¬ 
ticles and carry away a significant part of the collision 
energy. Thus, the effective energy for the production 
of secondary particles is the energy of interaction of a 
single quark pair, i.e. 1/3 of the entire nucleon energy. 
On the contrary, in the head-on heavy-ion collisions, the 
participating nucleons are considered colliding with all 
three constituent quarks from each nucleon. This makes 
the whole energy of the colliding nucleons (participants) 
available for the secondary particle production. Within 
this picture, one expects the results for bulk observables 
from head-on heavy-ion collisions at the c.m. energy per 
nucleon, to be similar to those from the pp/pp 

collisions but corresponding to a three times larger c.m. 
energy, i.e. at ^fSpp ~ 3^/snn- Such a universality is 
found to correctly predict 0 the value of the inidrapid- 
ity density in pp interactions measured at the TeV LHC 
energies . In addition, the multiplicity measurements 
in pp/pp interactions up to TeV energies are shown to be 
well reproduced by e“''e“ data as soon as the inelastic¬ 
ity is set to ~0.35 Q, i-e. effectively 1/3 of the hadronic 
interaction energy. This is in agreement with the dis¬ 
sipation energy picture where the structureless colliding 
leptons are considered to deposit their total energy into 
the Lorentz-contracted volume, similarly to nucleons in 
head-on nuclear collisions [2|. For recent discussion on 
the universality of hadroproduction up to LHC energies, 
see 0. 

Combining the above-discussed ingredients of the con¬ 
stituent quark picture and Landau hydrodynamics, one 
obtains the relationship between charged particle rapid¬ 
ity density per participant pair, p(ri) = {2/Npg,^t)dNch/dr] 
at midrapidity (rj fa 0), in heavy-ion collisions and in 
pp/pp collisions: 


p(0) _ 2iVch / Lpp 

Ppp{0) ~ fVpart NfP V ’ 


tSpp — 3-^/Saw ■ 


( 1 ) 


In Eq. (m , the relation of the pseudorapidity density and 
the mean multiplicity is applied in its Gaussian form as 
obtained in Landau hydrodynamics. The factor L is de¬ 
fined as L = \n{y/s/2m). According to the approach 
considered, m is the proton mass, Wp, in nucleus-nucleus 
collisions and the constituent quark mass in pp/pp colli¬ 
sions set to ^rup. Nch and are the mean multiplici¬ 
ties in nucleus-nucleus and nucleon-nucleon collisions, re¬ 
spectively, and A^part is the number of participants. 

Solving Eq. o for the multiplicity A^ch at a given ra¬ 
pidity density p(0) at -xJsnn, and for the rapidity density 
Ppp{0) and the multiplicity N^// at 3^/snn, one finds: 


^-^ch _ Tupp p(0) L 2 In 3 

A^part Ppp(O) y ln(4.5y'SArAr/mp) ’ 

y/sNN = y/s/f/,/3 . (2) 

Further development, as outlined below, treats this de¬ 
pendence in terms of centrality [l|. The centrality is re¬ 
garded as the degree of the overlap of the volumes of the 
two colliding nuclei, characterized by the impact param¬ 
eter. The most central collisions correspond, therefore, 
to the lowest centrality while the larger centrality to the 
more peripheral collisions. The centrality is closely re¬ 
lated to the number of nucleon participants determined 
using Monte Carlo Glauber calculations. Hence, the 
largest number of participants contribute to the most 
central heavy-ion collisions. The centrality is thus related 
to the amount of energy released in the collisions, i.e. to 
the effective energy, Enn. The latter, in the framework 
of the proposed approach, can be defined as a fraction of 
the c.m. energy available in a collision according to the 
centrality, a: 


£nn = '/snn{1 — a). (3) 

Conventionally, the data are divided into centrality in¬ 
tervals, so that a is the average centrality per centrality 
interval, e.g. a = 0.25 for the centrality interval of 20%- 
30% centrality. 

Then, for the effective c.m. energy Enn, Eq- © reads 


2iVch ,,pp pjO) L 21n3 

A^part Ppp(O) y In {4.5 ENN/rrip) ’ 

£NN = y/s/f/,/3 , (4) 

where p{0) is the midrapidity density in central nucleus- 
nucleus collisions measured at ^/sj/// = Eatat. 

In fact, each of the scalings described by Eqs. m 
and (U) regulates a particular physics ingredient used 
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in the modelling of the PDE approach. The scaling 
introduced by Eq.([2|) embeds the constituent quark 
model, which leads to establishing a similarity between 
hadronic and nuclear collisions. The scaling driven by 
Eq.dS]) addresses the energy budget effectively retained 
in the most central collisions while considering the global 
variables from noncentral collisions. 


III. MULTIPLICITY C.M. ENERGY 

DEPENDENCE IN CENTRAL NUCLEAR AND 
pp/pp COLLISIONS 

Pigure[T] shows the c.m. energy dependence of the mul¬ 
tiplicity measured in head-on nucleus-nucleus collisions 
(solid symbols) in the energy range of ^snn = 2 GeV 
to 2.76 TeV. Given the fact that the measurements sup¬ 
port the second-order logarithmic dependence on y/sJm 
up to the top RHIC energy @,[il while the power-law 
dependence is obtained for the LHC data , we fit the 
head-on data by the “hybrid” fit function, 

2 N 

= (-0.577 ± 0.177) -k (0.394 ± 0.094) ln(sArAr/so) 

fVpart 

-k(0.213± 0.014) ln^(sArAr/so) 

-f (0.005 ± 0.009) (sArAr/so)(° ’^®=^°-^^^(5) 

Here, so = 1 GeV^. In the following, the factor so is taken 
the same in all fit functions and omitted for brevity. This 
fit is shown in Pig. [T] by the solid line. Note that from 
the theoretical description point of view, the logarithmic 
dependence is considered to characterize the fragmenta¬ 
tion source(s) while the power-law behavior is believed to 
come from the gluon-gluon interactions [32| ; for a review, 
see [13 ■ 

We also fit the head-on collision multiplicities with 
the power-law function. The power-law dependence of 
the multi plicity is expected in different theoretical ap¬ 
proaches |34l436l | and the data from nuclear and pp/pp 
collisions seem to follow this type of behavior at higher 
energies 0 , 113 . The power-law fit gives: 

= (-6.72 ± 1.44) + (5.42 ± 1.11) (6) 

-/Vpart 

and is shown by the dashed-dotted line in Fig. [TJ 

Recently, it was shown [131 that the multiplicity of the 
gluon-gluon interactions are better described within a 
nonequilibrium statistical relativistic diffusion model us¬ 
ing log^(sArAr) dependence. Using this behavior we fit the 
c.m. energy dependence by the corresponding fit func¬ 
tion, 

= (0.72 ± 1.85) -k (0.75 ± 0.39) ln(s7V7v) 

^ part 

-k(0.019 ± 0.002) In^(sArAr). (7) 


Here, the linear-log term reflects the multiplicity from 
the fragmentation sources, as noticed above. The fit is 
shown by the pink solid line in Fig. [T] The fit is made 
starting the lowest NA49 energy in order to match the 
LHC multiplicity. One can see that the fit seems to be 
indistinguishable from the power-law function for the en¬ 
tire fit range, and is slightly below the power-law behav¬ 
ior above the current LHC data. Some enhancement in 
the low-energy range is expected due to no gluonic source 
considered to be present at these energies. 

In addition to these fits, we show the log^(sArAr)-fit 0, 
01 up to the top RHIC energy (thin dashed line). One 
can see that the power-law fit well describes the data 
and is almost indistinguishable from the hybrid fit up 
to the LHC data. Some minor deviation between the 
two fits can be seen in the range from the top RHIC 
energy to the LHC energy. Meantime, the second-order 
log polynomial lies below the data for ^snn > 200 GeV. 
This observation supports a possible transition to a new 
regime in heavy-ion collisions at ^Jsnn of about 1 TeV, as 
indicated earlier in the studies of pseudorapidit]^article 
and transverse energy densities at midrapidity [l[. 

Addressing now Eq. (I2|), we calculate the mean mul¬ 
tiplicity A'ch/(Vpart/2) for nucleus-nucleus interactions 
using the pp/pp measurements. The calculated values 
are shown in Fig. [T] by large open symbols. The ra¬ 
pidity density PppiO) and the multiplicity N//^ are taken 
from the existing data Q or, where not available, calcu¬ 
lated using the corresponding experimental ^jSpp fits^ at 
y/Spp = 3 y/SNN, in accordance with the approach con¬ 
sidered here. The p(0) values are as well taken from 
the measurements in central heavy-ion collisions wher¬ 
ever available, while for the nonexisting data the experi¬ 
mental ht^ is used. 

One can see that the calculated Vch/(0.5Apart) val¬ 
ues follow the measurements from nucleus-nucleus col¬ 
lisions at y/sJm from a few GeV up to the TeV LHC 
energy. The observed agreement between the heavy-ion 
measurements of Ach/(Apart/2) and the values obtained 
from the pp-based calculations points to the universality 
of the multiparticle production process in different types 
of collisions. 

Solving Eq. for the mean multiplicity N/!// in pp 
collisions, we estimate its values for > 2 TeV to 

be about 47 at = 2.36 TeV, 67 at 7 TeV, and 79 
at 13 TeV with 5% uncertainties. Here for the calcula¬ 
tions, one uses the fit to the heavy-iou midrapidity den¬ 
sity data p(0), as described above, and the fit by ALICE 


^ The powewr-law fit, Eq. is used for , while Ppp{0) is 
calculated using the linear-log fit po^O) = -0.308-1-0.276 In(spp) 
[3 and the power-law fit by CMS [Sal . ppp(O) = —0.402-j-, 
at yjSpp < 53 GeV and at y/Spp > 53 GeV, respectively. 

^ The linear-log fit p(0) = —0.33 -j- 0.38 In(siVAf) 0,113 is applied 
at y/SNN ^ 63 GeV, and the power-law fit p(0) = 0.73 
applied above 63 GeV as recently reported by ALICE using the 
measurements up to = 5.02 TeV [33l . 
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FIG. 1; The energy dependence of the charged particle mean multiplicity per participant pair. The large solid symbols show 
the measurements from the most central nucleus-nucleus (AA) collisions given as a function of the nucleon-nucleon c.m. energy, 
^sjvAT. The calculations by Eq. © based on pp/pp data at the c.m. energy ^Spp = 3 ^/snn are shown vs ^Spp/3 by large 
open symbols. The small open symbols show the AA data at different centralities as a function the effective energy enn 
(Eq. ©). The RHIC centrality data are shown after removing the energy-balanced limiting fragmentation scaling ingredient, 
while this ingredient is not taken into account for the LHC centrality data (see text). The multiplicity data of the most-central 
AA collisions are measured by the ALICE experiment at LHC [l3|, by the PHOBOS experiment at RHIC BEi, by the 
NA49 experiment at CERN SPS [2^ and by the E895 experiment at AGS [2j| (for the latter see also [l^). The low-energy 
HADES measurements at GSI are taken from [^ . The centrality data are taken from the measurements by the PHOBOS 
experiment at RHIC 0 and by the ALICE experiment at the LHC [ 13 , The values obtained from Eq. m for the AA 
mean multiplicity are based on: nonsingle diffractive pp collisions at FNAL by the E735 experiment [3, 0, at CERN by the 
UA5 experiment at ^yspp = 546 GeV and -^/Spp = 200 and 900 GeV [1^; pp collisions from CERN-ISR a nd from the 
inelastic data from the bubble chamber experiments [2^ - [3^ . the latter having been compiled and analyzed in j^. The LHC 
multiplicities in pp interactions are calculated using the hybrid fit obtained here, Eq. ®. The solid and the dashed-dotted show, 
correspondingly, the hybrid fit, —0.577-f 0.394 In(sjviv) + 0.213 In^(sjviv) + 0.005s^|f^, and the power-law fit, —6.72-|-5.42 
to the most central AA data. The thin dashed fine shows the second-order log fit —0.35 -I- 0.24 In(sivjv) -|- 0.24 In^(sjviv) to the 
most central AA data up to the top RHIC energy UM- The dashed and the dotted fines show, correspondingly, the hybrid fit, 
3.04— 1.41n(ejviv) + 1.121n^(ejvjv) -f 0.032and the power-law fit, —6.62 + 5.43to the centrality AA data. The pink 
solid line and the thin dashed-dotted line show the fits 0.72+0.75 In(sivjv)+0.019 ln®(sjvjv) and 1.7+2.36 ln(£jviv)+0.152 ln®(£jvjv) 
to the most central collision and centrality A A data, respectively (see text). The right-inclined hatched area shows the prediction 
for heavy-ion collisions at ^sjviv = 5.13 TeV and the left-inclined hatched area gives the prediction expected from pp collisions 
at ^Spp — 13 TeV. 
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FIG. 2: The c.m. energy dependence of the charged particle mean multiplicity in pp/pp collisions. The measurements are taken 
from: nonsingle diffractive pp collisions at FNAL by the E735 experiment [tI. [^. at CERN by the UA5 experiment at y/Spp = 
546 GeV [5^ and y/Spp = 200 and 900 GeV [2^; pp collisions from CERN-ISR [27||. and from the inelastic data from the bubble 
chamber experiments [2^4311 .1 ^ . The solid line shows the hybrid fit, 1.60—0.03 ln(spp)+0.18 ln^(spp)+0.03 s^'p^, the dotted line 
shows the power-law fit, —7.36 + 6.97 and the dashed line shows the second-order log fit, 3.18 — 0.57 In(spp)+ 0.216 In^(spp). 

The open stars show the predictions by the PDE approach with the error bars of 5% uncertainty. 


to the head-on heavy-ion data on the mean multiplicity 
m [similar to the results for the fits of Eqs . ([S]) and dH])] , 
along with the LHC measurements (dlMT^ of the pseu¬ 
dorapidity density Ppp(O). The calculated values of 
are shown as a function of ^Spp by open stars in Fig. O 
along with the existing multiplicity measurements from 
pp/pp collisions. 

The measured dependence on ^/Spp in the en¬ 
ergy range spanning the interval between a few GeV to 
1.8 TeV are fitted with the power-law, second-order log- 
polynomial and the hybrid functions. The hybrid and 
the power-law fits read, 

= (1.60 ± 0.23) + (-0.03 ± 0.10) In(spp) 

+ (0.18 ± 0.01) ln2(spp) + (0.03 ± 0.02) (8) 

and 

Atpp = (-7.36±0.16) + (6.97±0.12)4p^^^=^°°°^\ (9) 

correspondingly. 

From Fig. [5] one can conclude that the available data 
do not give any preference to one or another fit func¬ 
tion. This is similar to the pre-LHC observations where 


the power-law fit were found to be indistinguishable from 
the log^ polynomial fit at > 53 GeV Q. Interest¬ 
ingly, these two functions are also found to fit equally 
well the nonsingle diffractive midrapidity density, as ob¬ 
tained by CMS: cf. fits in [d^ vs those in [s^. The fit 
functions start to deviate from each other at the c.m. en¬ 
ergy above a few TeV but still not far one from another 
even at ^JSpp ^ 10 TeV. This may point to apparently 
no change in the multihadron production in pp interac¬ 
tions up to the highest LHC energy, in contrast to a new 
regime possibly emerging at •s/snn ~ 1 TeV in heavy-ion 
collisions. 


It is remarkable how well the PDE predictions on 
-^ch y/^pp > 2 TeV follow the power-law fit made 
to the measurements at y/Spp <1.8 TeV. This and the 
above-indicated “no change” in the hadroproduction in 
pp collisions as soon as one moves to TeV energies are in 
an agreement with the prediction 0 which seems to be 
the only successful one for the midrapidity density in pp 
collisions at = 7 TeV [T^ . 
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FIG. 3: The charged particle mean multiplicity per participant pair as a function of the number of participants, Apart- The 
solid circles show the dependence measured in AuAu collisions at RHIC by the PHOBOS experiment at ^/snn = 62.4, 130 
and 200 GeV [l^ (bottom to top). The solid stars show the measurements from PbPb collisions at the LHC by the ALICE 
experiment at y^sTcv = 2.76 TeV [2^. The triangles show the calculations by Eq. o using pp/pp data. The dashed lines 
represent the calculations within the effective-energy approach based on the hybrid fit, Eq. m, to the c.m. energy dependence 
of the mean multiplicity in the most central heavy-ion collisions shown in Fig. [T] The dashed-dotted line show the predictions 
for ^/sNN = 5.02 TeV using the average of the fits, Eq. IS) and Eq. ®. The open squares show the effective-energy calculations 
which include the energy-balanced limiting fragmentation scaling (see text); the solid lines connect the calculations to guide 
the eye. The open circles show the PHOBOS measurements at ^/snn = 200 GeV multiplied by 2.87. The open stars show the 
ALICE measurements at -^snn = 2.76 TeV multiplied by 1.3. 


IV. MULTIPLICITY CENTRALITY 
DEPENDENCE 


In this section, we address the point whether the cen¬ 
trality dependence of the mean multiplicity from heavy- 
ion experiments is described by Eq. similarly to the 
midrapidity pseudorapidity density in [l|. In Fig. [31 we 
show the Apart-dependence of Ach/(Apart/2). The data 


are taken from the measurements by the PHOBOS ex¬ 
periment at RHIC and by the ALICE experiment at 
LHC [13 . The solid triangles show the estimations using 
Eq. dH). As above, in the case of the ,/sATAr-dependence, 
the rapidity densities Ppp(O) and p(0), and the multiplic¬ 
ity are taken from the existing data (3 or, where not 
available, are calculated from the fits described above. 
According to the consideration developed here, p(0) is 
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taken at -^/snn = £nn, and Ppp{0) and are taken at 

^/^pp — 3 ^ NN • 

One can see that the calculations, which are driven by 
the centrality-defined effective c.m. energy snn, well re¬ 
produce the LHC data except slightly underestimating 
a couple of the most peripheral measurements. For the 
RHIC data, however, the difference between the calcula¬ 
tions and the measurements is visible already for medium 
centralities, i.e. for more central collisions. These ob¬ 
servations are also interrelated with the difference ob¬ 
served in the measurements at RHIC vs those from LHC. 
Indeed, at RHIC, the participant-pair-normalised mean 
multiplicity is found to be independent of centrality, 
while a decrease with centrality, or monotonic increase 
with A^part, is observed at the LHC. This becomes even 
clearer when the 200 GeV PHOBOS data are multiplied 
by a factor of 2.87 (open circles in Fig. [3]) which allows 
matching the ALICE data from the highly central colli¬ 
sions. 

In Fig. 131 the above-obtained c.m. energy fit, Eq. m, 
made to the head-on collision data, is applied to the cen¬ 
trality measurements at -^snn = e nn and the results are 
shown by the dashed lines.^ The observations made for 
the calculations are valid here as well. This points to 
the complementarity of central collisions and centrality 
data once the calculations are made in the c.m. effective- 
energy EjvAT terms. 

To clarify the observed differences, in the following 
sections the distributions of the pseudorapidity density 
are investigated in the context of the PDE picture 
considered here. 


V. PSEUDORAPIDITY DENSITY 
DISTRIBUTION IN HEAVY-ION COLLISIONS: 
CENTRAL VS NONCENTRAL COLLISIONS 


Figure [T] shows the distributions of charged particle 
pseudorapidity density per pair of participants measured 
in head-on and very central heavy-ion collisions and in 
pp/pp interactions. The heavy-ion data represent the 
PHOBOS measurements made in AuAu collisions at the 
RHIC at = 19.6, 62.4 and 200 GeV [Ii| and the 

ALICE measurements from PbPb collisions at the LHC 
at y'sAfAT = 2.76 TeV da. The distributions from pp/pp 
interactions are taken as measured by the UA5 experi¬ 
ment [i^ at .Js/p/ = 53 at the ISR and at = 200 
GeV at the SPS, by the P238 experiment at the SPS 


[46l| and by the CDF experiment at the Tevatron m 


= 630 GeV, and by the CMS [43 , LHCb 


48| 


and TOTEM experiments at 


= 7 TeV at the 


LHC. The data shown are taken at the c.m. energies 
Spp ~ 3 ^ysNN or Senn- Except for the CMS mea- 


^ Similar results are obtained from Eq. ©. 


surements, the negative -?7 data points from pp/pp inter¬ 
actions are the reflections of the measurements taken in 
the positive-/; region. 

Within the considered model of constituent quarks 
and the Gaussian form of the pseudorapidity distribu¬ 
tion in Landau hydrodynamics, the relationship between 
the pseudorapidity density distributions p(ri) and Ppp(r/) 
reads 


p(??) _ 2 Vch ^ 21n3 _- 77 ^_ 

Pppiv) fVpartA^fh V Lnn ^ Latat (2-l-LArAr/ ln3) 

( 10 ) 

Here, all variables are defined the same way as in Eq. o, 
i.e. taking into account the constituent quark scaling of 
the c.m. energy as soon as one relates pp/pp interactions 
to central heavy-ion collisions. 

Using Eq. m, the heavy-ion distributions are calcu¬ 
lated based on the Ppp{rf) spectra shown in Fig. |31 The 
calculated distributions are shown by solid symbols in 

Fig.m 

One can see that the calculations are in very good 
agreement with the measurements. Minor deviations are 
due to some mismatch between ^Spp and (or 

^£nn) and, as expected, due to a slight noncentrality; 
this is especially visible at y/SNN = 19.6 GeV where the 
energy mismatch is of a largest fraction. It is notice¬ 
able how well the PDE picture allows one to reproduce 
the pseudorapidity density distributions from heavy-ion 
interactions in the full -77 range, from central -77 to forward- 
77 regions, in the ^^snn range spanning over more than 
2 orders of magnitude. Remarkably, the pseudorapid¬ 
ity density distributions, measured in pp/pp collisions, 
despite being above those measured in heavy-ion colli¬ 
sions at 19.6 GeV or, on the contrary, lying far below the 
heavy-ion data from the LHC almost in the full -77 range, 
equally well reproduce the heavy-ion data as soon as be¬ 
ing recalculated within the PDE approach. Interestingly, 
the calculations at the LHC energies, well reproduce the 
heavy-ion data despite the pp measurements from the 
three different experiments are combined. A slight devi¬ 
ation in the negative -77 region is due to some asymmetry 
in the ALICE data. 

Let us now address peripheral collisions to clarify the 
deviation in centrality dependence between the data and 
the calculations as it is observed in Fig[3l 

In Fig. [51(a), the distribution p{r]) measured 0 in 
AuAu collisions by the PHOBOS experiment at i/satat = 
130 GeV at 45a = 0.475, is shown along with the Ppp{0) 
distribution measured in pp collisions by the UA5 exper¬ 
iment at y/Spp = 200 GeV [i^, i.e. at ^Spp Senn 
according to our approach. 

Applying Eq. (ITUl) . we calculate the p{r]) spectrum 
which is shown in Fig. jSKa) by solid squares. The calcu¬ 
lations agree well with the measurements in the central- 
77 region while fall below the data outside this region. 
This finding shows that in noncentral collisions, the cal¬ 
culations within the approach, which combines the con- 
























FIG. 4: The pseudorapidity distributions of charged particle pseudorapidity density per participant pair. The open circles 
show the measurements by the PHOBOS experiment in AuAu collisions at RHIC at (a) -^snn = 19.6 GeV, (b) 62.4 GeV 
and (c) 200 GeV p^ . and (d) by the ALICE experiment in PbPb collisions at the LHC at ^/snn = 2.76 TeV [^. The open 
triangles show the distributions measured in pp interactions by the UA5 experiment at ^/Spp = 53 GeV at the ISR and at 
^yspp — 200 GeV at the SPS [1^, by the P238 experiment at the SPS and by the CDF experiment at the Tevatron at 
^yspp = 630 GeV, and in pp interactions by the CMS [d^l, LHCb [4^ and TOTEM experiments at -^Spp = 7 TeV at the 
LHC. The solid markers show the calculations by Eq. using pp/pp data at ^Spp « 3 ^/snn or Ssnn- Apart from the CMS 
data, the negative-r; data points for pp/pp interactions are the reflections of the measurements taken in the positive-?; region. 


stituent quark picture and the relativistic Landau hy¬ 
drodynamics, reproduces well the pseudorapidity density 
around the midrapidity while underestimate the mean 
multiplicity. The former conclusion is well confirmed by 
our recent studies reported in [l[ for the midrapidity ob¬ 


servables, and the latter one is demonstrated by Fig. [H 
To clarify the obtained features, the following com¬ 
ments are due. 

In the PDF picture proposed here, the global observ¬ 
ables are defined by the energy of the participating con- 
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FIG. 5: (a) The charged particle pseudorapidity density per participant pair as a function of pseudorapidity. The open circles 
show the distribution measured in AuAu collisions at RHIC by the PHOBOS experiment at y/SNN = 130 GeV in 45%-50% 
centrality interval, fl^ . The open triangles show the distributions measured in pp interactions by the UA5 experiment at 
the SPS at yJSpp — 200 GeV [45|| . The solid squares show the distribution calculated from Eq. GSI by using the UA5 pp 
data at ^Spp ~ Senn (see Eq. ([3| for the definition of enn)- The solid circles show the beyond-midrapidity part obtained 
from the calculations using the energy-balanced limiting fragmentation scaling, i.e. under the shift r; —>■ 77 — ln{E nn / s nn)■ 
The negative -77 data points for pp interactions are the reflections of the measurements taken in the positive -77 region, (b) 
Same as (a) but the measured distributions of AuAu and pp collisions are shifted by the beam rapidity, p' = rj — T/beam, with 
ybeam = \n{^/s/mp), where s is, correspondingly, snn or Spp, and the calculated distribution is shifted to 77 ' = 77 — pes with 
7 /eff = ln(ejviv/ 7 np). The distribution measured in AuAu collisions and the calculated distribution coincide in the fragmentation 
region, when being shifted by j/beam for AuAu data and by j/efi for the calculations, that represents the energy-balanced limiting 
fragmentation scaling. 


stituent quarks pumped into the overlapped zone of the 
colliding nuclei. Hence, the bulk production is driven by 
the initial energy deposited at zero time at rapidity 77 = 0 , 
similar to the Landau hydrodynamics. Then, as is ex¬ 
pected and mentioned above, the pseudorapidity density 
at midrapidity is well reproduced for all types of nuclear 
collisions, from the most central to peripheral ones. As 
shown in [ij, similarly, the centrality dependence of the 
transverse energy density at midrapidity is well repro¬ 
duced by the calculations and complements the c.m. en¬ 
ergy dependence of the head-on data. Note that this sim¬ 
ilarity in the pseudorapidity density and the transverse 
energy pseudorapidity density is in accordance with the 
same functional form of the (pseudo)rapidity density dis¬ 
tribution obtained either in the longitudinally expanding 
system considered in the original Landau model or when 
the development in the transverse direction is included 


From Fig. [SKa), one can see that the calculated 
distribution p{r]) is narrower than that of the data. The 
narrowness of the calculated distribution with respect 
to the measured one is explained by a smaller value of 
Snn compared to the value of the actual collision energy 


However, the calculations in Eq. (flUl) are made 
with the multiplicity Nch taken from the most central 
nucleus-nucleus collisions at the c.m. energy equal to 
Enn (and similarly in Eqs. ([2]) and (|3]) for the midra¬ 
pidity density p(0)). In other words, in the approach 
applied here, similar to the Landau hydrodynamics, the 
collisions of nuclei are treated head-on-like. 


VI. ENERGY-BALANCED LIMITING 
FRAGMENTATION 

It is established that at high enough energies, in differ¬ 
ent types of interactions the pseudorapidity density spec¬ 
tra, measured at different c.m. energies, become similar 
in the fragmentation region. It means that they are inde¬ 
pendent of a projectile state (which is the beam or target 
rest frame) for the same type of colliding objects, i.e. be¬ 
ing considered as a function of 77 ' = r/ — where 

2 /beam = ln(y^s)^/ 77 T.p) is the beam rapidity [a, [l3]- This 
observation obeys a hypothesis of the limiting fragmen¬ 
tation scaling |53l |. 
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Considering the limiting fragmentation hypothesis 
within the effective-energy approach, one expects the 
limiting fragmentation scaling of the distribution /o(? 7 ), 
which is measured at to be similar to that of the 

calculated distribution but taken at the effective energy 
Enn- Note that the limiting fragmentation phenomenon, 
though being expected as an universal phenomenon for 
the Gaussian form of p(r7) [l3,[13,[ll|: naturally arises in 
Landau hydrodynamics [9|. 

In Fig. [Sljb), the limiting fragmentation hypothesis is 
applied to both the measured and the calculated pseudo¬ 
rapidity density distributions p{r]) from Fig. [5ja) using 
the c.m. energy and the effective energy, respectively. 
Therefore the measured distribution p(r]) is shifted by 
the beam rapidity, j/beam, while the calculated distribu¬ 
tion from Eq. (HOD is shifted by peS = hi(eNN/'mp) and 
becomes a function of rj' = rj—yes, as expected. One can 
see that the calculated p{r]') distribution of noncentral 
heavy-ion collisions agrees well with the measured distri¬ 
bution p{r/'). This finding points to a new energy scaling 
as soon as the effective-energy approach is applied. In 
analogy with the limiting fragmentation scaling, we call 
the observed scaling the ^^energy-balanced limiting frag¬ 
mentation scaling”. Due to this scaling, the calculated 
pseudorapidity density is getting corrected outside the 
central -?7 region accordingly. 

To this end, in Fig. Elja), the calculated distribution 
p{r]) is shifted by the difference (j/eff — ybeam) in this re¬ 
gion: 77 —>■ ?7 — nn / y/s nn) , or, using the effective 

energy definition, Eq. ([3]), p ^ p — ln(I — a). The cal¬ 
culated distribution p{p)^ where the shift is applied, is 
shown by the solid circles in Fig. [SJa). The shift bal¬ 
ances the energy and this brings the calculations to the 
measured pseudorapidity density distribution in the full- 
p range in noncentral heavy-ion collision. It is clear that 
in head-on or very central collisions, a approaches zero 
which makes the shift negligible (cf. Fig. 2]). 

This finding allows obtaining Nch within the PDE ap¬ 
proach. Namely, the difference between the two Nch val¬ 
ues, one obtained by integrating the calculated pseudo¬ 
rapidity density distribution from Eq. (nni, and another 
one of the same distribution but being shifted to the left 
by ln(I — a), is added to the iVch value obtained from 
Eq. (|4]). Where no pseudorapidity density distributions 
are available in pp/pp measurements at ^/Spp = Seat/v, 
the energy-balanced limiting fragmentation scaling is ap¬ 
plied to reproduce the calculated p(r]): the measured dis¬ 
tribution from a noncentral heavy-ion collision is shifted 
by (j/beam “ Ves), i-s. p ^ p -\- ln(l - a). Then 
is calculated as above, by adding to the calculation of 
Eq. Q the difference between the integral from the ob¬ 
tained shifted distribution and the measured multiplicity 
in this noncentral heavy-ion collision. 

Using this ansatz, the values of Wh are calculated for 
each centrality for the RHIC measurements. The cal¬ 
culations are shown by open squares in Fig. [3] One can 
see that now the calculations well reproduce the measure¬ 
ments from RHIC, with no deficit in noncentral collisions. 


The energy-balanced limiting fragmentation scaling 
provides an explanation of the “puzzle” between the cen¬ 
trality independence of the iVpart-normalized mean mul¬ 
tiplicity and the monotonic decrease of the normalized 
midrapidity pseudorapidity density with the centrality, 
as observed at RHIC. As shown above, the pseudorapid¬ 
ity density at midrapidity is determined by the effective 
energy of centrally colliding nucleon participants. Hence, 
the value of this observable increases towards head-on 
collisions as soon as the effective energy, made available 
for particle production, increases with increasing number 
of participants (decreasing centrality). However, the mul¬ 
tiplicity is measured in the full 77 -region, so it gets addi¬ 
tional contribution from beyond the midrapidity. In the 
context of the PDE picture, this contribution is due to 
the balance between the collision c.m. energy shared by 
all nucleons of colliding nuclei and the centrality-defined 
effective energy of the interacting participants. The more 
peripheral is the collision, the larger the additional con¬ 
tribution is. This contribution can be directly estimated 
by the energy-balanced limiting fragmentation scaling, 
introduced here, which leads to scaling between the mea¬ 
sured pseudorapidity distribution and distribution calcu¬ 
lated within the PDE approach. 

From Fig. [31 one can conclude that, in contrast to the 
RHIC measurements, almost no additional contribution 
is needed for the PDE calculations of Eq. (|4]) in order to 
describe the LHC mean multiplicity data. As the calcu¬ 
lations imply, they are made by considering the nucleus- 
nucleus collisions as head-on collisions at the c.m. energy 
of the value of Enn (p( 0) in Eq. (|4]) as well as Nch in 
Eq. (TTOll are taken from the head-on ^/snn fits). How¬ 
ever, as shown above, the additional contribution to the 
mean multiplicity increases with increasing collision cen¬ 
trality, i.e. while going towards more peripheral collisions. 
For head-on collisions, however, this contribution tends 
to zero. Given the multiplicity measurements at the 
LHC are well reproduced without the energy-balanced 
additional contribution, one concludes that in heavy-ions 
collisions at the LHC at TeV energies the multihadron 
production obeys a head-on collision regime, for all the 
centrality intervals measured. This points to apparently 
different regimes of hadroproduction occurring in heavy- 
ion collisions with ^/sJm between a few hundred GeV 
and TeV energies. This observation supports a similar 
conclusion made above, which is suggested from the ob¬ 
servation of a change of the functional type of the fit 
needed to describe the energy behavior as soon as the 
LHC data are included, see Fig. |TJ 

The discussed difference between the mean multiplic¬ 
ity, and hence the full pseudorapidity density distribu¬ 
tions, measured in heavy-ion collisions at RHIC and at 
LHC have been earlier addressed in [13 , where the 

model of three sources, the gluon-gluon midrapidity and 
two quark-gluon fragmentation sources, are applied to 
understand the observations from experimental data. In 
the context of the PDE approach given here, the differ¬ 
ence in the nature of collisions at effective c.m. energy is 
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appealed to explain different centrality dependence of the 
data from the two colliders. Similarly to calculations in 
[Hil, additional contribution from the fragmentation 
regions are shown to be needed at RHIC. However, no 
such contribution is needed at the LHC energy. Mean¬ 
time, the midrapidity pseudorapidity densities measured 
at RHIC and at LHC do not show different behavior with 
centrality and are found [I| to be similarly well repro¬ 
duced by the PDE calculations where no preference is 
given to midrapidity or fragmentation sources. 

There are other approaches, which also consider 
the three effective regions in pseudorapidity density 
distributions of charged particles produced in pp/pp 
and in heavy-ion collisions. In the string percolation 
model , the fragmentation region is populated by 

strings of valence quarks and the midrapidity region by 
additional short strings between quarks and antiquarks. 
In other approaches, one introduces a leading particle 
activity within the hydrodynamic [57l - [6?il | or thermal [6l| 
pictures of the multiparticle production processes. Like 
in the PDE consideration, presented here, a similarity 
of the mechanism of particle production in pp/pp and 
heavy-ion collisions is also assumed in these approaches. 
Whereas within the PDE approach the leading particles 
resulting from the spectators are considered to be pro¬ 
duced in nucleon-nucleon collisions, where a single quark 
pair interaction is assumed, no leading particle effect 
is implied for central nucleus-nucleus collisions, where 
the entire energy of the participants is considered to be 
available for bulk hadron production. As already noticed 
above, also no difference in the particle production 
sources in different pseudorapidity regions is assumed 
in the PDE approach. Then, the c.m. energy scaling 
due to the key picture of the constituent quarks, applied 
to the Landau hydrodynamics, allows revealing the 
universality of the multihadron dynamics in hadronic 
and nuclear interactions. 


VII. MULTIPLICITY EFFECTIVE ENERGY 
DEPENDENCE IN HEAD-ON AND 
NONCENTRAL COLLISIONS 


Given the obtained agreement between the data and 
the calculations, and considering the similarity put for¬ 
ward for £ MM and ■^smn, one would expect the measured 
centrality data at Emn to follow the y/SMN dependence of 
the mean multiplicity in the most central nuclear colli¬ 
sions. In Fig. [U the measurements of the charged particle 
mean multiplicity of head-on nuclear collisions are added 
by the centr ality measurements by the PHOBOS [l^ and 
the ALICE experiments (Fig. [3]) where the central¬ 
ity data are plotted as a function of Emn- Due to the 
above finding of the energy-balanced limiting fragmenta¬ 
tion scaling, explaining the lack of centrality dependence 
of the mean multiplicity at RHIC energies, these data are 
plotted by subtracting the energy-balanced contribution. 


In addition, the centrality data at -^Jsmn = 19-6 GeV are 
shown in Fig. [1] but not in Fig. [31 From Fig. [I] one con¬ 
cludes that effective-energy dependence of the centrality 
data complements the c.m. energy behavior of the head- 
on collision data. 

To better trace the similarity between the head-on col¬ 
lision and centrality data, we fit the £ArAr-dependence of 
the centrality data by the hybrid and the power-law func¬ 
tions, similarly to the head-on collisions. For the hybrid 
fit one gets 


2 N 

= (3.04 ± 0.60) - (1.40 ± 0.24) Iiy^Emn) 

-^part 

+{1.12 ± 0.04) In^(eiVAr) -f (0.032 ± 0.028)e^mn'^^^°-^°%1) 

This fit is shown by the dashed line in Fig. [TJ This fit 
agrees well with the same type of the fit to the head- 
on collision data in the entire available energy range 
though lying slightly above the latter one for the data 
at y/SNN ~ 10 GeV. For the log^(£ArAr) fit function of the 
three-sources approach, similar to Eq. o, one finds: 

2 iV 

= (1.70 ± 1.49) -f (1.18 ± 0.54) ln(£ArAr) 

^part 

-f (0.152 ± 0.008) In^iEMN)- (12) 

The fit is shown by the thin dashed-dotted line in Fig. lU 
and lies on top of the analogous fit, Eq. (O to the head-on 
data, except a slight enhancement at y/SMN ~ 10 GeV, 
similar to the hybrid fit. The power-law £ 7 V 7 V-fit for the 
centrality data is shown by the dotted line in Fig.jTJ It is 
found to be similar to the power-law ^SAw-fit, Eq. ([6|), 
to the head-on collision data shown by the dashed-dotted 
line. 

From this one concludes that within the picture pro¬ 
posed here the data are well reproduced under the as¬ 
sumption of the effective energy which governs the mul¬ 
tiparticle production. This points to the the same en¬ 
ergy behavior in multihadron production for all types of 
heavy-ion collisions, from peripheral to the most central 
collisions. 

Here, let us stress an important corollary of the PDE 
approach. As soon as the effective energy in nucleus- 
nucleus collisions determines the pseudorapidity density 
at midrapidity, then the midrapidity pseudorapidity den¬ 
sities at the same effective energy but at different c.m. 
energy get the same value. In other words, the densities 
are defined by the effective energy independent of the en¬ 
ergy of the collision. The observation made here for the 
multiplicity dependence on the effective energy confirms 
the observation made earlier for the midrapidity densi¬ 
ties [l| , while it adds another important ingredient which 
takes into account the additional energy-balanced con¬ 
tribution to the mean multiplicity in noncentral nucleus- 
nucleus collisions. 

From the hybrid fits obtained, we estimate the mul¬ 
tiplicity for the future LHC heavy-ion run. Since the 
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hybrid fit for the head-on collision data and the fit to 
the centrality data show slightly different increase with 
c.m. energy, the predictions of the two fits are averaged. 
Hence, the mean multiplicity 2iVch/iVpart value is pre¬ 
dicted to be about 119 with 5% uncertainty in the most 
central heavy-ion collisions at y/SNN = 5.13 TeV. The 
prediction is shown by the right-inclined hatched area 
in Fig. [TJ This value is close to the value of about 116 
one gets from the ALICE fit El and about 108 which 
one obtains from Eq. ([7]). In addition, the fit-averaged 
prediction based on pp collisions at ^Spp = 13 TeV, re¬ 
calculated within the PDE approach, is shown in Fig. [1] 
as the left-inclined hatched area. 

The predictions are made as well for the centrality de¬ 
pendence and are given in Fig. [31 We give the predic¬ 
tions for = 5.02 TeV, considering the recent mea¬ 

surements reported by ALICE for the rapidity density 
and expecting the mean multiplicity measurements 
at this energy. The two types of predictions are shown. 

First, similar to the above predictions made to the 
head-on collisions, we use the fit functions. As soon as, 
within the effective-energy approach, we treat noncen¬ 
tral collisions as central collisions at energy Enn, then 
we use the head-on collision multiplicity fits to predict 
the centrality dependence. This is similar to that made 
for the existing data as shown in Fig. [3] by the dashed 
lines. However, for the predictions, we use the average 
values of the hybrid and the power-law fits, Eqs. ([S]) and 
as soon as those deviate for ^/snn above 2.76 TeV. 
The prediction for ^^snn = 5.13 TeV centrality depen¬ 
dence is shown by the dashed-dotted line. The centrality 
and Vpart values are alike in the 2.76 TeV data shown. 
The expectations show an increase of the mean multi¬ 
plicity with iVpart (decrease with centrality) from about 
52 to about 118. The increase looks to be slightly faster 
than at ^/snn = 2.76 TeV, especially for the peripheral 
region. 

Second, the PDE set of prediction is made using the 
calculations based on Eq. (|3|) combining the constituent 
quark model and the Landau hydrodynamics. This pre¬ 
diction for the centrality dependence is shown by the solid 
inverted triangles. One can see that the predictions are 
close to ones obtained from the head-on collision data 
fits. 

The predictions are compared with the LHC data at 
y/SNN = 2.76 TeV. To better match the predictions for 
highly central collisions, the 2.76 TeV data points are 
multiplied by 1.3. One can see that the predictions are 
well reproduced by the scaled data. This indicates no 
change of the hadroproduction mechanism expected with 
increase of the c.m. energy at LHC, in contrast to what 
is seen as one compares the scaled 200 CeV RHIC data 
with the 2.76 TeV measurements. 

An interesting issue to be addressed in the framework 
of the PDE picture is asymmetric collisions, such as 
nucleon-nucleus (p/d-nucleus) ones. In these interactions 
the multiplicity is also expected to have no centrality 
dependence. This is due to the many nucleon-nucleon 


interactions of the incident proton with the nucleons 
of the interacting nucleus while the secondary particles 
produced in the reaction are assumed to be created 
out of the c.m. energy deposited to the interaction 
zone. The proton and nucleus are considered to interact 
via a single pair of constituent quarks, one from the 
proton and another one from a nucleon in the inter¬ 
acting nucleus. Then, no centrality dependence of the 
multiplicity is expected in p-nucleus collisions with the 
multiplicity values to be similar to that from pp/pp 
interactions at the c.m. energy ^yspp ~ As a 

consequence, this, at a given centrality, results in Vpart/2 
for the ratio Nch/N^^. These features have been indeed 
obtained in d-Au interactions at RHIC [6^. Moreover, 
the effect of the Vpart-dependence of the multiplicity 
ratio obtained at RHIC has been also observed in 
hadron-nucleus collisions at lower y/sJm ~ 10-20 GeV 
[^ . These observations seem to be also obtained 
at LHC, where the dependence of the pseudorapid¬ 
ity density on ,ysNN measured in p-Pb collisions at 
= 5.02 TeV is observed to be consistent with 
a power-law ^Spp-&t to pplv p data and its centrality 
dependence is shown [65], l66l | to demonstrate the im¬ 
portance of nucleon-nucleon interactions for p-Pb results. 


VIII. SUMMARY AND CONCLUSIONS 

In summary, the multihadron production process in 
nucleus-nucleus collisions and its universality in nuclear 
and hadronic interactions are studied. The study exploits 
the charged particle mean multiplicity dependencies on 
the c.m. collision energy per nucleon, yJSNN: and on 
the number of nucleon participants, or centrality, mea¬ 
sured in the energy range of a few GeV to a few TeV. 
The study is carried out in the framework of the earlier 
proposed approach of the dissipating effective energy of 
constituent quark participants ii, or the participant 
dissipating energy (PDE) approach. In this approach, 
the participants are considered to form the initial zone of 
a collision and to determine the production of hadrons at 
the very early stage of the collision. In this consideration 
one combines the constituent quark picture with Landau 
hydrodynamics and interrelates the multihadron produc¬ 
tion in different types of collisions by a proper scaling of 
the c.m. energy of collisions. In particular, an energy¬ 
scaling factor of 1/3 in pp/pp measurements is shown to 
reveal the universality of the multiplicity dependencies in 
nucleon-nucleon and nucleus-nucleus interactions. 

In the entire available range of about a few TeV, 

the energy dependence of the multiplicity in head-on col¬ 
lisions is found to be well described by the calculations 
performed within the effective-energy approach based on 
pp/pp data. Meanwhile, depending on the data sample, 
the calculations are found either to describe the mea¬ 
sured centrality dependence or to show some deviation 
between the calculations and the data. For the RHIC 
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data, the deficit in the predictions is observed for non¬ 
central collisions so that the predictions do not follow a 
constancy with the centrality as it is observed at RHIC. 
The LHC mean multiplicity centrality dependence, how¬ 
ever, is found to be well described by the calculations 
including the increase towards the most central collisions. 

To clarify the observations, approach of the effective- 
energy of the quark participants is applied to the pseudo¬ 
rapidity density distribution measured in heavy-ion colli¬ 
sions. The energy-balanced limiting fragmentation scal¬ 
ing is introduced based on assumption of the similarity 
of the fragmentation region of the measured distribu¬ 
tion in the beam rest frame and that determined from 
the calculations by using the effective energy. The re¬ 
vealed scaling allows us to reproduce the pseudorapid¬ 
ity density distributions independently of the centrality 
of collisions and then to correctly describe the central¬ 
ity independence of the mean multiplicity measured at 
RHIC. Moreover, this finding provides a solution to the 
RHIC “puzzle” of the difference between the centrality 
independence of the mean multiplicity vs the monotonic 
decrease of the midrapidity pseudorapidity density with 
the increase of centrality. The mean multiplicity is shown 
to get a fraction of additional contribution to account for 
the balance between the collision c.m. energy shared by 
all nucleons and the effective energy of the participants. 
However,the midrapidity pseudorapidity density is fully 
defined by the effective energy of colliding participants. 

Given that the calculations made in the context of the 
proposed approach are considering central collisions of 
nuclei, an agreement between the calculations and the 
LHC data indicates that at TeV energies the collisions 
seem to present head-on collisions of the participants at 
the c.m. energy of the scale of the effective energy. Thus, 
no energy-balanced additional contribution is required 
even with a relatively small number of participants at 
TeV energies. 

Based on the above findings, the complementarity of 
the head-on collisions and the centrality data is shown 
resulting in the similar energy behavior of the mean mul¬ 
tiplicity measurements as soon as the data are considered 
in terms of the effective energy. A departure of the c.m. 
energy dependence of the data from the second-order log¬ 
arithmic behavior to the power-law or higher-order loga¬ 
rithmic polynomial function one observes at the region of 
1 TeV suggests a transition to a new regime in nucleus- 
nucleus collisions at TeV energies. Interestingly, these 


findings made for a full collision rapidity range are similar 
to those drawn from the studies [l| of the pseudorapidity 
density and the transverse energy density at midrapidity. 
This is also in accordance with the change of the multi¬ 
plicity dependence on centrality which also indicates a 
possible change of the regime of multihadron production 
as one moves from the RHIC to LHC energies. 

The hybrid and the power-law fits are found to describe 
well the existing data on the charged particle multiplicity 
from pp/pp interactions in the entire c.m. energy ^Spp 
range up to the top Tevatron energy of 1.8 TeV. However, 
in this case no clear change from the power-law behavior 
to the quadratic log polynomial one is obtained in the 
multiplicity c.m. energy dependence. Moreover, the pre¬ 
dictions made here for the mean multiplicity for ^Spp in 
the LHC energy range of 2.36 to 13 TeV within the PDE 
approach demonstrate a closeness between the predicted 
values and the lower-energy yfSpp fit. One concludes that, 
in contrast to heavy ions, no change in multihadron pro¬ 
duction in pp collisions is expected up to the foreseen 
LHC energy. 

Based on the results of the hybrid fits, the predictions 
for the charged particle mean multiplicity in head-on 
heavy-ion collisions at ^/snn = 5.13 TeV at the LHC are 
given. Within the obtained complementarity of head-on 
collisions and centrality data, the predictions are made 
for the mean multiplicity centrality dependence to be 
measured. 

The soon-to-come measurements at the LHC are of 
crucial importance for further understanding of the mul- 
tihadron dynamics. This will shed the light on the univer¬ 
sality of the multihadron production process in different 
types of collisions and clarify the PDE approach and the 
obtained energy-balanced limiting fragmentation which 
have been shown to successfully describe the features of 
global key observables by relating hadronic and nuclear 
collisions. 
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